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Abstract. A simplified transient energy-transport system for semiconductors subject 
to mixed Diriclrlet-Neumann boundary conditions is analyzed. The model is formally 
derived from the non-isothermal hydrodynamic equations in a particular vanishing mo- 
mentum relaxation limit. It consists of a drift-diffusion-type equation for the electron 
density, involving temperature gradients, a nonlinear heat equation for the electron tem- 
perature, and the Poisson equation for the electric potential. The giobal-in-time existence 
of bounded weak solutions is proved. The proof is based on the Stampacchia truncation 
method and a careful use of the temperature equation. Under some regularity assump- 
tions on the gradients of the variables, the uniqueness of solutions is shown. Finally, 
numerical simulations for a ballistic diode in one space dimension illustrate the behavior 
of the solutions. 



1. Introduction 

The basic model for the charge transport in semiconductor devices are the drift-diffusion 
equations for the electron density and the electric potential. This model gives fast and 
satisfactory simulation results for devices on the micrometer scale, but it is not able to cope 
with so-called hot-electron effects in nanoscale devices. A possible solution is to incorporate 
the mean energy in the model equations, which leads to energy-transport equations, first 
presented by Stratton [IB] and later derived from the semiconductor Boltzmann equation 
by Ben Abdallah and Degond [1]. The analysis of the energy-transport model is very 
involved due to the strong coupling and temperature gradients. Therefore, we consider in 
this paper a simplified energy-transport model which still includes temperature gradients 
but the coupling to the energy equation is weaker than in the full model. An important 
feature of our model is that it is derived formally from the hydrodynamic semiconductor 
equations in a zero relaxation time limit, which provides a physical modeling basis without 
heuristics (see Section [2]). Our goal is to prove the existence and uniqueness of solutions 
to this model and to provide some numerical illustrations of the solutions. 
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The model consists of a drift-diffusion- type equation for the electron density n{x,t), a 
nonlinear heat equation for the electron temperature 6{x,t), and the Poisson equation for 
the electric potential V{x,t): 

(1) dtn-div{\/{ne) +nVV) =0, 

77 

(2) dw{K{n)ve) = -{e-eL{x)), 

T 

(3) -X'^AV = n- C{x) inn, t> 0. 

Here, ^(n) is the thermal conductivity, 9l{x) the lattice temperature, and C{x) the doping 
profile characterizing the device under consideration. The scaled physical parameters are 
the energy relaxation time r > and the Debye length A > 0. Equations ([I])-© hold in 
the bounded domain C M'^ (rf > 1) with the initial condition 

(4) n{0) = ni in Vt. 

We suppose that the boundary dVt G C°'^ consists of two parts To and Fat satisfying 
dVL = r/)UrAr, ri:inrAr = 0, FatIs closed, and the {d — l)-dimensional Lebesgue measure 
of To is positive. The electron density, temperature, and potential are assumed to be known 
on the Dirichlet boundary, which models the contacts, whereas the Neumann boundary 
models insulated boundary parts: 

n = ud, = Od, V = Vd on Td, 

Vn-u = Ve-u = VV-u = Q onTN, 

where v denotes the exterior unit normal vector on dVL. 

Before we detail our analytical results, we review related models in the literature. First, 
temperature effects have been included in the drift- diffusion equations by allowing for 
temperature-dependent diffusivities [T7j or temperature-dependent mobilities [HI [121 [201 [23] 
coupled to a heat equation. Typically, the so-called non-isothermal drift-diffusion equations 
are of the form 

dtn — div Jn = 0, J„, = DVn + /iraVV, 

dtO - dw{K{e)ve) = F, F = j^-vv + w, 

where J„ is the particle current density, D and fi are the diffusivity and mobility, respec- 
tively, and W = —n{9 — Ol)/t is the relaxation term. The difficulty in these models is 
that the Joule heating term J„ ■ W involves quadratic gradients of the potential, which 
resembles the thermistor problem; see, e.g., [21]. However, temperature gradients in 
which need to be taken into account, have been ignored. 

In Xu allowed for temperature gradients in J„ but he truncated, as in [23], the 
Joule heating term by setting F = max{0, J„ ■ W + W} in order to allow for a maximum 
principle. A different approach was adopted in [3] , where a kind of quasi- Fermi potential 
via = nexp{—V/9) was introduced. Ths model of [3] includes temperature gradients, 
but the coefficient contains the electric potential which is not the case in the energy- 
transport models derived in [1]. We also mention non-isothermal systems with simplified 
thermodynamic forces which were studied in [1]. 
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Compared to our model ([I])-®, the energy-transport equations contain cross- diffusion 
terms also in the energy equation A typical form of these models reads as 

dtn-dwJn = 0, j;„ = V(nr) 

^dtinO) - divJe = Jn-VV + W, Je = V(n^"+^) + n^"VV, 

where the parameter a > is related to the elastic scattering rate in the collision operator 
(see Example 6.8 in [I5]). In our model (II])-(I2]), a = 1, and the diffusion scaling implies 
that the variation of the energy density, ^dt{n6), and the Joule heating term are negligible 
(see Section [2]). The main difficulty of the above model is that the corresponding diffusion 
matrix is neither diagonal nor tridiagonal and that it degenerates for = or ^ = 0. 
Existence results were achieved for stationary equations near thermal equilibrium [TUl [H] 
and for the transient model El El] if the initial data are close to the stationary drift- 
diffusion solutions. General existence results, both for the stationary and time-dependent 
model, were proved in [TJ [S] but the diffusion matrix was assumed to be uniformly positive 
definite, thus avoiding the degeneracy. All these results give only partial answers to the 
well-posedness of the problem, and a complete global existence theory for the energy- 
transport equations for any data and with physical transport coefficients is still missing. 

In this paper, we wish to bring forward the existence theory for energy-transport-type 
models by analyzing the system ([I])-([3]), whose complexity is in between the well-understood 
drift-diffusion model and the energy-transport equations. In fact, in our model, the energy 
equation simplifies such that the application of the maximum principle for 6 becomes pos- 
sible. The remaining difficulties are due to the drift term nVO in ([T]) and the quasilinearity 
K{n) in ([2]). Note that, in view of the mixed boundary conditions, we cannot expect the 
regularity G L°° which would simplify the existence proofs significantly. 

Our main idea is a careful use of the temperature equation in order to deal with the 
drift term nV6. More precisely, we replace this term in ([1]) formally by 

div(nV^) = div (-nVe) = -div(KVe) +Vn-Ve- -VO ■ ( ^Vn + 

and using (|2]), we find that ([1]) equals 

(6) dtn - div(^Vn) = div(nVV) + —{9 - 9l) Vn-V9- -^\V0\\ 

The computations will be made rigorous on a weak formulation level in Section [3l From 
the above formulation we see that the last term on the right-hand side models a sink if 
di%/d9 > 0. This condition is satisfied, for instance, in the case of the Wiedemann- Franz 
model K{n, 9) = n9. By the maximum principle, we expect to obtain an upper bound for 
n. 

However, we need the stronger condition On/ 89 = 0. The reason is that the lack of time 
regularity for 9 makes it difficult to deal with nonlinear terms, such as 6'Vn, to prove the 
continuity of the fixed-point operator. Although in physical models, it is often assumed 
that the thermal conductivity depends on the temperature 9, a dependency on n only 
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also occurs in the physical literature. For instance, the choice K{n) = n was suggested in 
[T^ Formula (2.16)] to study spurious velocity overshoots in hydrodynamic semiconductor 
models. 

From the physical application, we expect that the electron density n stays positive if 
it is positive initially and on the Dirichlet boundary parts. Even if k depends on n only, 
the proof of a positive lower bound for n is not obvious, since it is not clear how to deal 
with the term Vn ■ in ([6]) which is in only. We suppose that either K{n) is strictly 
positive or n{n) = n. In the former case, we avoid any degeneracy; in the latter case, 
{n/ n){dn/ dn) = 1, and the term involving Vn • V6 in (j6]) vanishes. 

Motivated by the above considerations, we impose the following conditions on the ther- 
mal conductivity: Let k, G (^""^([O, oo)) such that there exist Kq, Ki, n^,, n* > with 

(i) k{z) > for all 2; > 0, 
(7) (ii) either k{z) > kq > for all z > 0, or k{z) = z for all < 2; < n^=; 

(iii) K,{z) > HiZ for all z > n*. 

Condition (i) allows for the degenerate case k{0) = 0. Condition (ii) ensures the uniform 
ellipticity of equation ([2]). Indeed, if ^(n) = n for n < n*, we are able to prove that the 
solution n is strictly positive and then, ^(n) is strictly positive, too. The last condition is 
needed to prove an upper bound for the particle density. 
The boundary data are assumed to satisfy 

no, VDeL\0,T;H\n)), 60 e L^{0,T;W'''^{n)), 
no, Od G L°°(0,T;L°°(fi)), infri^ > 0, inf^z? > 0, 

where q > 2 and Qj- = Q x [0,T). The initial data and the given functions fulfill the 
conditions 

(9) m, Ol, C e L°°{n), infn/>0, inf 0^ > 0, inf C(x) > 0. 

S~2 S~2 S~2 

In order to deal with the mixed Dirichlet- Neumann conditions, we introduce the space 

H^{n U r^) = {ue H\n) : « = O on Yd}. 
For properties of this space, we refer to (THl Chapter 1.7.2]. Furthermore, we set H^^{Q U 

Theorem 1 (Existence of solutions). Let Q G M.'^ (d > 1) be a bounded domain with 
dQ G C^'^, T, T, \ > 0, and let k G C"'^([0, 00)) satisfy ([7]). Furthermore, assume that 
([HD and Q hold. Then there exists a weak solution in,e,V) G L'^{0,T; H^{fl)f to (HD-dSD 
satisfying dfU G L^(0,T; H~^{Q U F^v)) and 

< n{t) < Koe^\ < m < e{t) <M mVt, t^ (0, T). 

Furthermore, if k{z) = z for < z < n^, 

n(t) > A;oe-°* > inQ,te{0,T). 
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In the above theorem, the constants are defined by 

JiTo = max < n*, supn/, sup ud, supC(x) 
I n rDx(o,T) n 



(10) 



ko = min < n^, inf n/, inf riu 

I n rox(0,T) 



M = max<sup0L, sup 0£)>, m = min < inf inf Od 

n rox(o,r) J I ^ rox(o,T) 

1 1 « ^ 

a = -SUpt^L + "TT, P 



The proof of the theorem is based on the Leray-Schauder fixed-point theorem and the 
Stampacchia truncation method. In particular, the truncation is needed in the diffusion 
coefficients of div(6'Vn) and div(fi;(n)V^) to make these expressions uniformly elliptic. 

Due to the quasilinearity of the temperature equation, we are able to show the uniqueness 
of solutions only in a function space which includes bounded gradients. 

Theorem 2 (Uniqueness of solutions) . Let the assumptions of Theorem [I] hold and let 
K he locally Lipschitz continuous on [0,oo). Then there exists a unique solution {n,6,V) 
to ([I])-(l5]) in the class of bounded weak solutions satisfying n G L°°{Q,T]W^'^{VL)), 9 G 
L°°(0,T;Pyi'~(fi)), where p > 2 if d = 2 and p > d if d > 3. 

The paper is organized as follows. Equations ([I])-(l3]) are formally derived from the 
hydrodynamic model in Section [2l The existence theorem is proved in Section |3l and 
Section H] is devoted to the proof of the uniqueness theorem. In Section |5l we present 
numerical results for a simple one-dimensional ballistic diode illustrating the behavior of 
the electron temperature in the presence of a cooling and heating lattice temperature. 

2. Derivation of the model equations 

Equations ([I])- ([3]) are formally derived from the (scaled) hydrodynamic model (see, e.g., 
[151 Chapter 9]): 

dfU — div Jn = 0, 

5^ j„ _ div ( III^] - V(n^) - nVV = 

\ n J Tp 

n f 3 

dt{ne) - div( J„(e + d)) - J„ ■ W - div(fi;(n, 9)V9) = e - -^^ 

Te \ 2 

and V is given by the Poisson equation (j3]). Here, J„ denotes the particle current density, 
is a tensor product, Tp is the momentum relaxation time, and Te the energy 
relaxation time. The energy density is the sum of the thermal and kinetic energies: 

ne = -nO H 

2 2n 
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Energy-transport equations can be derived from the vanishing momentum relaxation 
limit. To this end, we set e = Tp and rescale the equations by t — )■ t/e and J — )■ eJ. This 
corresponds to the physical situation of a long time scale and small current densities. The 
rescaled equations become: 

3 |J„P 

(11) dm — div .L = 0, ne = -nO H 

(12) e'dtJr, - e^diY (^ '^"^'^" ^ - V{ne) - nVV = - J„, 

n / 3 

edt{ne) - £div( J„(e + 9)) - eJn ■ VV - div(fi:(n, d)\/d) = i e - -9l 

re \ 2 

In the formal limit e — 0, we obtain the limiting model 

77 / 3 \ 3 
dtn-dw{V{ne) +nVV), diY{K{n,9)V9) = - ie - -9l] , e = -9, 

which corresponds to ([I])- ([2]) with r = 2re/3. 

In the literature, usually a different limit is performed in order to derive energy-transport 
equations. Indeed, if we rescale additionally k ^ en (small thermal conductivity), and 
assume that the energy relaxation time is of the same order as the momentum relaxation 
time, r = tq = e, the rescaled energy equation reads as 

(13) edt{ne) - ediv( J„(e + 9)) - eJ^ ■ VV - ediv(/t(n, ^) V^) = ~en (^e - ^9 

Then, dividing this equation by e and performing the formal limit £ — > in ( !TT|) and (fT2ll . 
we find the usual energy-transport model with particular diffusion coefficients (see p!5l 
Chapter 6.4]). 

Our simplified model is valid in diffusive situations in which the thermal conductivity is 
strong and the energy relaxation time is much larger than the momentum relaxation time. 

3. Proof of Theorem [T] 

The existence proof is based on the Leray-Schauder fixed-point theorem and a truncation 
method. For this, we consider the truncated problem 

(14) dtn - div(ft„,MVn) = div(n^V(^ + V)), 

(15) div{^{nk,K)V9) = '^{9-9L), 

T 

(16) -X'^AV = nK- C{x) inVt, t> 0, 
with the initial and boundary conditions (|l])-([5]), where 

nK = max |0, minjK, n}}, 
nk,K = max | fc, mm{K, n} } , 
9m,M = max |m, min{M, 9} } , 
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and k = k{t) = fcoe~"*, K = K{t) = Kqc^^. We recall that the constants m, M, ko, Kq, a, 
and (3 are defined in fITU]) . Observe that the lower truncation of n in f[T51) is not necessary 
if K{n) > Ko > for all n > 0. In this case, we replace K,{nk^K) by k^uk)- 
We divide the proof in several steps. 

Step 1: Definition of the fixed-point operator. Let w G L^(0, T; and a G [0,1]. 

For given t G (0,T), let V^(t) G be the unique solution to the linear problem 

-X'^AV{t) = w{t)K - C{x) in fi, V{t) = Voit) on T^, VV^(t) ■ = on T^v- 

Since w G L^(0, T; L^(f2)), we find that V : (0,T) — )■ H^{Q) is Bochner-measurable and 
1^ G L\0,T-H\n)) (see, e.g., [2l pp. 1133f.]). 

Next, let 6'(t) G H^{Q) be the unique solution to the linear uniformly elliptic problem 

div {K{w{t) k,K)^e) = ^^^^^{6 -Ol) in 6 = 9D{t) on To, ■ z/ = on T^. 

Again, the integrability of w allows us to conclude that 6 G L^(0,T; H^{Q)). 
Finally, consider the linear parabolic problem 

dtn - div(6'„_MVn) = a-div(wxV(6' + V)) in fi, t> 0, 

n = ano on To, Vn ■ u = on Tn, n{0) = anj in Q. 

Since the right-hand side of the parabolic equation is an element of L^(0,T; H^^IQUTm)), 
there exists a unique solution n G L^(0,T; H^{fl)) fl H^{0,T; H^^^Q U Fat)). This shows 
that the operator S : ^^(0, T; L^iQ)) x [0, 1] ^ ^^(0, T; ^^(l])), (w, or) h-> n, is well defined. 
It holds that S{w, 0) = for all w G L2(0, T; 

By using 6' — 6*/) as a test function in (fT5|) . standard estimates and the lower bound of k 
show that 

||^l|L2(o,r;Hi(Q)) < Ci, 

where Ci > depends on kq, m, M, K, 9l, and 9d- Similarly, 

where C2 > depends on i^. A, C{x), and Vd- Therefore, employing n — an^ as a test 
function in ([T]), a Gronwall estimate implies that 

\\n\\L2(o,T;m(n)) + \\dtn\\L2(o,T;H-'^{nnr,^)) < C3, 

where C3 > depends on m, K, nn, nj, ci, and C2. 

We claim that ^ is slightly more regular. Indeed, using the (admissible) test function 
{6 - M)+ = max{M, 6} in ([TS]), we obtain 

[ |V(e-M) + pdx < / K(nfc,i^)|V(e- M)+|^da: 
Jn Jn 

= -- [ nKi9-9L){9-My <0, 
T Jn 

since 9 — 9l > on {9 > M}, where = > or = min^gj^ hi{z) > (see ([71)). We 
infer that ^ < M on f2, t > 0. In a similar way, the test function {9 — m)~ = min{m, 9} 
yields 9 > m. In particular, we have 9m,M = 9. Thus, the right-hand side of the heat 
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equation is an element of L°°(0,T; L°°{Q)). By elliptic regularity, we have [131 Theorem 1] 
9{t) e W^'P{n) for some 2 < p < g, and hence, G Lp(0,T; LP(fi)). 

Step 2: Continuity of the fixed-point operator. Let Wj — )■ w strongly in L^(0, T; L'^{fl)) 
and (Tj — )■ cr as j — )■ cxD. Let 6*^ and Vj be the solutions to 

(17) divWK)„)V9,) = M£(9, -9J, -A^AV-.^Kk-CW in a 

with the corresponding boundary conditions. Then, by the above elliptic estimates, up to 
a subsequence, 

Oj e, Vj V weakly in L^(0, T; H^{n)). 

Since hi{{wj)k,K) i^iwk^K) strongly in U'{Q,T] U'iyt)) for any r < oo, we can pass to the 
limit in (fT71) to obtain 

diY{K{wkK)ye) = —{6 -Ol), -\^^V = wk-C{x) infi. 
r 

In view of the compact embedding H^iVt) Aubin's lemma shows that L^(0,T; 

n H\0,T] H'\Q U Fjv)) is compactly embedded into L^{0,T; L^{Q)). Thus, the 
above estimate for Uj proves that, again up to a subsequence, 

rij n strongly in ^^(0, T; 

We have to show that n = S{w, a). This is proved by passing to the limit in the parabolic 
equation satisfied by Uj. The problem is the limit of (OjVuj) since {6j) and (Vrij) both 
converge only weakly. We claim that 6j — )■ 6 strongly in L'^{0,T; H^{Q)). Taking the 
difference of the equations satisfied by 9j and 9, respectively, and using 9j — 9 as a. test 
function, we find that 

[ I K{{w,\K)\y{9,-9)\^dxdt 
Jo Jn 

= - I I {<{'^i)KK) - K{Wk,K))'^9 ■ V{9j - 9)dxdt 
Jo JQ 

- - / / - ^kW + K)a'(% -9)- iiw,)K - wk)9l) (9, - 9)dxdt 

T Jo Jn 

< - / / {<iwj)k,K) - n{wk,K))V9 ■ V{9.j - 9)dxdt 
Jo Jn 

- - / / {{wj)K - WK)9{9j - 9)dxdt 
T Jo Jn 

+ - I I - WK)dLi9j - 9)dxdt. 

Jo Jn 

The regularity G L^{0,T; L^{Q)) for some p > 2 and the strong convergence of 
ti{{wj)k,K) -> i^{wk,K) in any U{Q,T]U{Vt)) imply that {K{{wj)k,K) - ^■{wk,K))^9 
strongly in ^^(0, T; ^^(fi)). Hence, since V^^ V9 weakly in ^^(0, T; L2(fi)), the first 
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integral on the right-hand side converges to zero. Similarly, in view of the L°° bounds for 
9 and 9^, the second and third integrals converge to zero. Since hi{{wj)k,K) > > 0, this 
shows the claim. 

Hence, we can pass to the limit j — oo in the equation 

I {dtnj,(l))dt+ [ I 9jVnj -Vcpdxdt = -aj I [ {wj)K'^{9j + Vj) ■ V(j)dxdt, 
Jo Jo Jn Jo Jn 

where (■, ■) denotes the dual product on Hq{QUT^) and G L'^{0, T; HQ{QUTiy)), to infer 
that n solves 

dtn - div(^Vn) = -adiY{wKV{9 + V)) in ^^(0, T; H~\n U T^)). 

This implies that n = S{w,a). Hence, S is continuous and, by the Aubin lemma, also 
compact. 

We prove uniform estimates in L°°{0, T; L°°{Q)) for all fixed points of S{-, a) which allows 
to remove the truncation and which yields uniform estimates in L^(0, T; L^(f2)) needed for 
the fixed-point theorem. 

Step 3: L°° bounds for n. Let n be a fixed point of 5'(-,cr). First, observe that the 
test function n~ in f[T^ immediately implies that n~ = and n>0inr2, t>0, since 
rix'Vn' = inQ. To derive an upper bound, we set u = e~^^n. Then u solves the equation 

(18) dtu - div(^VM) = adiv(Mi^„V(^ + V)) - (3u, 

since Uk = max{0, minji^oc''*, e''*M}} = e'^* max{0, minjii'o, «}} =: ^^^Ukq- Let L > K 
and define = n{nk,K)~^UKo{uL — Ko)~^y where ul = min{L,M}. This truncation is 
necessary to obtain G L'^{0, T; Hq^QUTn))- Furthermore, 0(0) = since u{0) = rij < Kq 
in Vt. We employ the test function in the temperature equation (1151) : 

(19) - / K{nk,K)V9-V(t)dx = - [ -!^(9-9L)uK,iuL-Koydx. 

Jn r Jn f^{nk,K) 

First, we compute the left-hand side: 

- / K{nk,K)V9 ■ V(f)dx = - UK,y9 ■ V{ul - Ko^dx - / Vma'o ■ V9{ul - Ko)^dx 
Jn ' Jn Jn 

+ / -^^Au,-KoypVnK-V9dx. 
Jn i^[nk,K) on 

The second and third integrals vanish since Vmxq = and Vn^ = on {m > Kq\. We 
obtain 

- / K{nk,K)'^9 ■V(f)dx = -Ko / V9 ■ V{ul - Ko^dx. 
Jn ' Jn 

Therefore, since 9 < M and uk / t^ijik^K) < ^k^x/ i^i^k^K) < l/'^i (see ([7])), (IT^ becomes 

(20) -Ko f V9 ■ ViuL - KoYdx = - [ -^^{9 - 9l)uk,{ul - Ko^dx 

Jn ^ Jn i^{nk^K> 



M 
< 



/ uko{ul - Ko)^dx. 
Jn 
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Next, we use {ul — Kq)^ as an admissible test function in f[T8|) . An elementary compu- 
tation shows that 



F{s) = J\aL - KoVda > ^{{sl - i^o)+)' 



Therefore, since F{u{0)) = F{nj) = 0, 



1 



{dtU,{uL-Ko)+)ds= / {F{u{t))-F{u{0)))dx>- {{u{t)L - Ko)+) dx, 
'0 Jn ^ Jn 

where (■, ■) denotes the dual product on H\(^VM^ fq). This gives 

\ [ {iu{t)L-Ko)+Ydx+ [ [ e\ViuL-Ko)+\^dxdt 
^ Jn Jo Jn 

<-a / UKo'^ie + V) ■ V{ul - KoYdxdt - (3 I I u{ul - KQ^dxdt. 
Jo Jn Jo Jn 

By the Poisson equation ( fT6l) . 

- / uk.VV ■ V{ul - Ko)^dx = -Ko [ VV ■ V{ul - K^^dx 
Jn Jn 

= -X^'^Ko [ (riK - Cix))iuL - KoYdx < 0, 
Jn 

since u > Kq is equivalent to n > K and hence, rix — C{x) = K — C{x) > Kq — C{x) > 
on {u > Kq}, using the definition of Kq. Then, taking into account (!20|) . we find that 

I [ {{u{t)L-KQ)+Ydx + m [ [ \V{uL-KQ)+\^dxdt 
^ Jn Jo Jn 

/ / ^Koi^L — KQ)^dx — /3 / / u{ul — KQ)^dxdt 
Jo Jn Jo Jn 

/ / UKoiuL - KQ^dx - (5 I / UKoiuL- KQ^dxdt 
Jo Jn Jo Jn 

-P] f I UK,{uL-KQ)+dxdt = 0, 
/ Jo Jn 

by the definition of /3. We infer that [ul — Kq)~^ = for all L > Kq. Letting L — )■ oo, 
we obtain {u — Kq)~^ = and thus, n < K in Q, t > 0. As a consequence, uk = n, 
and any solution to ffT^ - flTBl) solves ([I])-®. Furthermore, the L°° bounds provide the 
uniform estimates needed to apply the Leray-Schauder fixed-point theorem. This proves 
the existence of solutions to ([I])-®. 

Step 4- Positive lower bound for n. Assume that n{z) = z for all < z < n*. We claim 
that under this condition, n possesses a positive lower bound. In view of the upper bound 



M 

< 

M 



< 



M 
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from Step 3, [n — k)~ , where k = /cqC""*, is an admissible test function in f fT^ yielding 

(21) i [ {n-kyitfdx + m [ [ \V{n - k)-\'^dxdt < -a [ [ nVO ■ V{n - k)~dxdt 
2 Jn Jo Jn Jo Jn 



cr [ [ nVV -Vin- k)~dxdt + a [ [ 
Jo Jn Jo Jn 



k{n — k) dxdt. 
We write the second integral on the right-hand side as 



-a {n- k)W ■ V(n - k^dxdt -a / kW ■ V(ra - k^dxdt 
Jo Jn Jo Jn 

= -- [ [ VV ■V{{n-k)-Ydxdt-a [ [ kVV ■ V {n - k)' dxdt 
2 Jo Jn Jo Jn 

= --^ / [{n- C{x)){{n - kyfdxdt - ^ / [ {n- C{x)){n - k)' dxdt 
2A^ Jo Jn ^ Jo Jn 

< T^WCh^in) [ [ {{n-k)-Ydxdt + ^ [ [ k[-{n-k)-]dxdt, 
Jo Jn ^ Jo Jn 

using the Poisson equation and n{n — k)~ < k{n — k)~ in Vt. 

In order to estimate the first integral on the right-hand side of (pTj) . we employ the test 
function {n — k)~ in (fT5|) . Then, since K,{n) = n for < n < k < ko < n^, 

- n{e - eL){n - k)~dx = - K{n)V6 ■ V{n - k)~dx = - nV6 ■ V(ra - k)~dx. 
Jn Jn Jn 

Therefore, (12T!) becomes 

^ ^ {n - kyit^dx + m [ [ \V{n - k)-\^dxdt < - [ [ n{e - Oyin - k)' dx 

Jo Jn T Jo Jn 

+ ^l|C|U-{n)^ y {{n - k)-f dxdt + - [ [ k[-{n - k)-]dxdt 



Jn 



<^l|C|U^(n) r / {{n-krYdxdt 



1„„ „ 1 



Jn 

t 



+ + - «J J J k[-{n-k) ]dxdt 

+ - j j k[-{n - ky]dxdt = 0. 
We obtain (n — k)~ = and hence, n > k in Q, t > 0. 

4. Proof of Theorem [2] 

Let (ni,0i,Vi), {n2,92,V2) be two solutions to ([I])-© with the regularity indicated in 
the theorem. 
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Step 1: Estimate of"V{6i — 62) ■ We employ the test function 6*1 — 6*2 in the difference of 
the weak formulations for 6*1, 6*2, respectively: 

(22) f [ n{n2)\V{ei - e2)\^dxdt = - [ [ {n{ni) - K(n2))V^i ■ V(^i - e2)dxdt 
Jo Jn Jo Jn 

- - / / - 02) + {ni - n2){ei - eL)){ei - e2)dxdt. 

Jo Jq 

Using the Cauchy-Schwarz, Poincare, and Young inequalities, the second integral is esti- 
mated from above by 

C\\ni — n2||L2(o,T;L2(n))||^l " ^2 ||l2{0,T;L2(q)) 

< £:||V(6'i - 6'2)||i2(o^T;L2(n)) + - '^2||i2(o_T;L2(Q)). 

where c{e) > depends on e, the L°° bounds for di and 6^, and the Poincare constant. 
The Lipschitz continuity of k on [0, K] implies that 

(^(ni) - K(n2)) V^i ■ V(^i - 92)dxdt 

< c 



Jn 

[ [ \ni-n2\\Vei\\V{ei-e2)\dxdt 
Jo Jn 

< C\\ni - n2||L2(0,T;L2(n))i|V(6'i - 6^2) || l2(0,T;L2(o)) 

< £:||V(6'i - 6'2)||i2(o,T;L2(n)) + C{e)\\ni - n2||i2(o,T;L2(n)), 

where c{e) > depends on e and the L°° norm of V6'i. Since K{n2) > k^, > for some 
> 0, we find from (122!) . for e < k*/4, that 

(23) ||V(6'i - 02)\\L^O,T;L^n)) < c{K^)\\ni - n2||L2(0,T;L2(n)). 

Step 2: Estimate ofrii —n2 ■ We employ rii — 77-2 in the difference of the equations satisfied 
by rii and n2, respectively: 

(24) l f (rii -n2)(t)'dx+ / / e2\V{ni - n2)\^dxdt 
2 Jn Jo Jn 

t 



[ f {61- ^2) Vni ■ V{ni - n2)dxdt 
Jo Jn 

I [ (n2V(^i - 62) + {rii - n2)Vei) ■ V(ni - n2)dxdt 
Jo Jn 

[ [ {n2V{Vi - V2) + {rii - n2)VVi) ■ V{ni - n2)dxdt. 
Jo Jn 
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Applying Holder's inequality with p > 2 as in the theorem and l/p+l/g + l/2 = 1 to the 
first integral, we estimate as follows: 

- [ / (6*1 - e2)Vni ■ V{ni - n2)dxdt 
Jo Jn 

< \\6i — ^2||L2(0,T;L'J{C))||Vr2i||Loo(o^T;LP{C))||V(?T,i — ^22) || L2(0,r;L2(n)) 

< c\\V{6i - 6'2)||L2(o,T;L2(Q))||V(ni - n2)\\L^0,T;L^n)) 

< C\\ni - 7^2||L2(0,r;L2(f7))||V(ni - n2)||L2(0,T;L2(n)) 

< £||V(ni -n2)||i2(o,T;L2(n)) +C(^^)II^1 -^2||i2(o,T;L2(n))- 

In the second step we have used the Sobolev embedding H^{Q) and the Poincare 

inquality, and the third step follows from (1251) . 

For the second integral in (^^, we obtain, using again (12^ . 

/ / (ri2V(^i - 62) + (ni - n2)V^i) ■ V(ni - n2)dxdt 
Jo Jn 

< £||V(ni - n2)||i2(o,T;L2(n)) + C(£)||V(^1 - ^2)||i2(o,T;L2(n)) 

+ c(e)||ni -n2||i2(o,r;L2(n)) 

< e||V(ni - ri2)||i2(o,T;L2(n)) + 4^)11^1 - n2\\lHo,T;LHn))- 
Finally, for the third integral in fl24l) . we estimate 

f f {n2V{Vi - V2) + {ni - n2)VVi) ■ V{ni - n2)dxdt 
Jo Jn 

< e\\V{ni - n2)||i2(o,T;L2(f7)) + c(e)l|V(V^l - V2)\\l2^o,T;LHn)) 

+ ]- f [ Wi ■ V{ni - n2fdxdt. 
2 Jo Jn 

By the elliptic estimate for the Poisson equation, 

/ / (ri2V(Vl - V2) + (ni - n2)VVi) ■ V(ni - n2)dxdt 
Jo Jn 

< e|| V(ni - n2)||i2(o,T;L2(o)) + c{e)\\ni - n2||i2(o,T;L2(Q)) 
{rii — C{x)){ni — n2Ydxdt 



1 



2A 



Jn 



< e\\V{ni - n2)||i2(o,T;L2(n)) + C{e)\\ni - n2||i2(o,T;L2(n))- 



Inserting these estimates in flMl) and observing that 62 is uniformly bounded from below, 
i.e. 62 > m > in Q, t > 0, we infer that 

^11(^1 -^2)(t)||i2(f^) +m||V(ni -n2)||i2(o,T;L2(n)) 

< 3e\\V{ni - n2)||i2(o,T;L2(Q)) + 4^)11^1 - ^2||i2(o,T;L2(f^)). 
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Then, choosing e < l/(3m), the Gronwall lemma allows us to conclude that {ni — n2){t) = 
in Q for t > 0. This proves the uniqueness of solutions. 

5. Numerical approximation 

In this section, we present numerical results for the simplified energy-transport model 
with K{n, 6) = n9 on the interval [0, 1]. The initial and boundary conditions are 

ni{x) = C{x) for x G fi, n(0, t) = C(0), n(l, t) = C(l), 

^(0, t) = Olio), 6(1, t) = ^^(1), V{0, t) = 0, 1/(1, t) = U for t > 0. 

Equations ([I])-© are discretized on an equidistant grid with grid points Xi = iAx, where 
Ax = 1/(A^ — 1). The time grid points are tk = kAt, where At > 0. We employ central 
finite differences in space and the trapezoidal rule in time. Then, with the approximations 
n^, 9^, and V^^ of n{xi,tk), 9{xi,tk), and V{xi,tk), respectively, the discretized equations 
become 



1 

2(Ax)2 
1 

2rAx)2 



Given {n^'^ , 9^~^ ,Vl^~^) ^ we find {n\,9^,V^) by solving the above nonlinear equations 
subject to the corresponding (Dirichlet) boundary conditions using Newton's method. 
We simulate a ballistic diode which is defined by the doping profile 

C(x) = 1 + 0.25(tanh(100x - 60) - tanh(100x - 40)), x G [0, 1]. 

The physical parameters are given in Table [H and the scaled quantities are defined by 



qCras.^L'^ ' ° ° m„ ' Y ksTo ' t* ' 

For the computations, we choose = 201 grid points and the time step size At = 1.25 x 
10-^ 

We wish to study the impact of different lattice temperatures. First, we choose a lattice 
temperature which is cooling the interior of the diode, i.e. 9l{x) = ^(x — |)^ + |. Figure [1] 
shows the electron density and electron temperature at various times for applied voltages 
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Parameter 


Value 


Physical meaning 


ks 


1.3807 X 10~^^kgm/s^K 


Boltzmann constant 


eo 


8.8542 X 10-^2 A^s^/kgm^ 


Vacuum permittivity 


mo 


9.11 X 10"^^ kg 


Electron mass at rest 


? 


1.602 X 10"^^ As 


Elementary charge 






Maximum doping concentration 


To 


300 K 


Device temperature 


L 


75 nm 


Device length 


rrin 


0.067 -mo 


Effective electron mass 




11.7 


Relative permittivity of GaAs 


To 


0.9 X 10-12 s 


Momentum relaxation time 




3.0 X 10-3 


Scaled squared Debye length 


r 


3.126 


Scaled energy relaxation time 


Ko 


4.88 X 10-2 


Heat transfer coefficient 



Table 1. Physical and scaled parameters. 



U = 0.2 V and U = 1.0 V, respectively. In both cases, the electron temperature converges 
to its nonhomogeneous stationary profile as t — oo. Since the profile is convex, equation 
02]) implies that the particle temperature is larger than the lattice temperature. The 
profile of the electron density follows the doping profile except for the large applied bias 
U = 1.0 V. In this situation, the electric force is sufficiently strong to deplete the charge 
carrier concentration close to the left boundary point. 

Figure [2] illustrates the behavior of the electron density and electron temperature when 
the lattice temperature is heating the diode, i.e. 9l{x) = | — 3(x — |)2. Again, the electron 
temperature converges to a nonhomogeneous steady state, and the behavior of the particle 
density is similar to the case of cooling temperatures. The current- voltage characteristic is 
very close to that one with constant temperature (not presented). It can be seen that only 
for very large voltages, the current density becomes slightly smaller due to an increasing 
thermal energy fraction. This shows that the influence of the temperature equation is not 
very important in a ballistic diode although there are significant temperature gradients. 
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Figure 2. Electron density and temperature in the ballistic diode with 
heating lattice temperature at voltages U = 0.2 V (top) and U = 1.0 V 
(bottom). 
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